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Abstract 

A theory is offered for a novel device, mesoscopic four-terminal SQUID. The stud- 
ied system consists of a mesoscopic four-terminal junction, one pair of terminals of 
which is incorporated in a superconducting ring and the other one is connected with 
a transport circuit. The nonlocal weak coupling between the terminals leads to ef- 
fects of phase dragging and magnetic flux transfer. The behaviour of a four-terminal 
SQUID, controlled by the external parameters, the applied magnetic flux and the 
transport current is investigated. The critical current and the current voltage char- 
acteristics as functions of magnetic flux are calculated. In the nonlocal mesoscopic 
case they depend not only on the magnitude of the applied flux but also on its sign, 
allowing measurement of the direction of the external magnetic field. 
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The four-terminal SQUID is a superconducting device based on a four-terminal 
Josephson junction [1]. The multiterminal Josephson junctions (MTJ) gener- 
alizes the usual (two-terminal) Josephson junctions [2] to the case of weak cou- 
pling between several massive superconducting banks (terminals). Compared 
with two-terminal junctions, such systems have additional degrees of freedom 
and the corresponding set of control parameters, preset transport currents 
and (or) applied magnetic fluxes. As a result, the current- or voltage-biased 
and the magnetic flux-driven regimes can be combined in one multiterminal 
microstructure. One of the realizations of a MTJ presents the set of super- 
conducting terminals, which are connected with the short dirty microbridges 
having a common centre (Fig. la). In such (conventional) MTJ the nonlinear 
coupling between individual bridges is due to the fact that the supercon- 
ducting order parameter in the common centre is a function of the currents 
through all the microbridges. See review of macroscopic quantum interference 
effects in conventional multiterminal microstructures in Ref . [3] . Another type 
of multiterminal Josephson junction is based on the weak coupling of bulk 
superconductors through the two-dimensional clean normal layer (Fig. lb). In 
such a mesoscopic ballistic 4-terminal junction [4] the nonlocal coupling of 
supercurrents is established due to the phase dependent local Andreev levels 
inside the weak link. The nonlocal weak coupling in mesoscopic MTJ leads to 
the effects of phase dragging and magnetic flux transfer [5,6]. In the present 
paper we consider the behaviour of a mesoscopic 4-terminal SQUID (Fig. 2) 
emphasizing the specific features of mesoscopic case. 
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The studied system (Fig. 2) consists of a mesoscopic four-terminal junction, one 
pair of terminals of which is incorporated in a superconducting ring and the 
other one is connected with a transport circuit. The system can be controlled 
by an external magnetic flux $ e and by the transport current I. Current J 
is the self-induced circulating current in the ring circuit with self-inductance 
L. The state of the 4-terminal junction is determined by the dynamical vari- 
ables, phases (pi of the complex off-diagonal potential A exp(i<^j) in the i-th 
terminal (i=l..A). In the 4-terminal case we have three independent phase 
differences. Phases 9 and are the phase differences between terminals 2-1 
and 3-4, respectively. Phase x — l/2(v?i + V2) — 1/2 (v?3 + ^4) is the phase 
difference between the ring and the transport current circuit. The phase is 
related to the observable quantity, total magnetic flux threading the ring $, 
4> = For given external parameters / and $ e , Gibbs potential for the 
4-terminal SQUID in terms of variables <p, 9 and x nas the form [6]: 



(A) _ 2e$e)2 

U((j>, 9, X ; I, $ e ) = n— % — - -19- cos(9) - k cos(0) 

n ( 6 <P n ■ . 0\ , , 

-2 I a cos — cos — p sin 2 sin 2" / COS ^' ^ ' 

The last three terms in Eq.(l) are the Josephson coupling energy. The coeffi- 
cients K,a,/3 in Eq.(l) determine the coupling between the different terminals. 
They depend on the geometry of the weak link and on the transparency D 
of S — N interfaces. In a symmetric case of square iV-layer and ideal trans- 
parency (-D=l) they equal to : k — 1, a — y/2 + 1, (3 = y/2 — 1. The current / is 
measured in units of I = 7r7 12 A (T) 2 /4eT c , where 712 is the inverse Sharvin 
resistance between the terminals 1 and 2. L = (2e/h)LI K is the dimensionless 
self-inductance. 
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The specific feature of a mesoscopic 4-terminal SQUID is the appearance in 
the potential U of the term, proportional to coefficient (3. Note, that even for 
a completely symmetrical mesoscopic MTJ we have (3 7^ 0, contrary to the 
conventional MTJ (Fig. la). In the latter case (3 identically equals zero for any 
symmetry of microbridge configuration. 

The minimization of U with respect to phases 9, 0, x gives the equations 



I = sin 9 + 



. e 4> e.j 
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cos x = sign 
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The system of Eqs.(2-4) describes the behaviour of a mesoscopic 4-terminal 
SQUID in the stationary state. It determines the magnetic flux in the ring 
$ = ^0 as function of applied magnetic flux $ e and transport current /.The 
nonlocal Josephson coupling {(3 ^ 0) leads to the new effects. As follows from 
Eqs.(2-4) the magnetic flux $ in the ring produces the phase difference 9 on 
current driven junction 1-2 at zero current / flowing from terminal 1 to ter- 
minal 2. Similarly the transport current I induces in the ring the magnetic 
flux $ in the absence of external flux <3> e . The influence of the phase difference 
between one pair of terminals on the phase difference between another pair of 
terminals is what we call phase dragging effect in the mesoscopic 4-terminal 
junction. As applied to mesoscopic 4-terminal SQUID this effect results in the 
dependence of the critical current I c (maximal value of / at which the system 
of equations (2-4) has the solution) not only on the absolute value but also 
on the sign of the applied magnetic flux <3> e . The steady states domain in the 
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(/, <3> e ) plane is shown in Fig. 3 for the case of small self-inductance L <C 1. 
The boundary of the domain , curve J c ($ e ) (solid lines in Fig. 3), is 2ir peri- 
odic, but due to the terms proportional to (5 in Eqs.(2-4), it is not invariant 
under the transformation $ e — > — $ e . The symmetry is restored if we simul- 
taneously change $ e on — $ e and I on —I. Note, that in conventional case 
(P = 0) the boundary of the steady state domain J c ($ 6 ) is symmetric with 
respect to the axes (J, <3> e ) (dashed line in Fig.3). Outside the steady state 
domain, the stationary solutions for (9, $) are absent and system goes to the 
nonstationary resistive regime. We have studied the dynamical behaviour of 
the 4-terminal SQUID in the frame of the heavily damped resistively shunted 
junction (RSJ) model [2]. The features of the dynamical behaviour of the meso- 
scopic 4-terminal SQUID are again affected by the terms proportional to (3, 
i.e. by nonlocal coupling. The current-voltage characteristics in the transport 
channel, V(I), (V is the time averaged voltage between terminals 2 and 1) are 
shown in Fig. 4 for different values of $ e . The voltage V(I) in applied magnetic 
flux $ e depends on the sign of the $ 6 as well as critical current I c . Thus, the 
mesoscopic 4-terminal SQUID can be used for direction-sensitive detection of 
weak magnetic fluxes. In addition, the system described by the potential (1) 
permits the existence of the bistable states even for negligible dimensionless 
self-inductance L — > 0, which can be exploited in the superconducting phase 
qubit design (see Ref. [7]). 
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Figure Captions 

Fig.l: The multiterminal Josephson junction. Bulk superconductors 1,..4 are con- 
nected with the crossed microbridges (a) or are weakly coupled through the 
rectangular of a two-dimensional electron gas (b). 

Fig. 2: The mesoscopic four-terminal SQUID. 

Fig. 3: The steady state domain for mesoscopic four-terminal SQUID in plane 

(I, $ e ) of the control parameters (solid line). $ e is measured in units h/2e. 

Dashed line corresponds to the conventional four-terminal SQUID. 
Fig. 4: The current-voltage characteristics of a mesoscopic 4-terminal SQUID for 

different values of applied flux $ e . $ = ^ is the superconducting flux 

quantum, L = 0. 
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Fig. 3. 
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